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Abstract 

We study a Hubbard model with a strong onsite impurity potential based on an improved fluctuation-exchange (FLEX) approx- 
imation, which we call the GV 7 -FLEX method. We find that (i) both local and staggered susceptibilities are strongly enhanced 
around the impurity. By this reason, (ii) the quasiparticle lifetime as well as the local density of states (DOS) are strongly sup- 
pressed in a wide area around the impurity (like a Swiss cheese hole), which causes the "huge residual resistivity" beyond the 
s-wave unitary scattering value. These results by the GV 1 method naturally explains the various impurity effects in HTSC's in a 
unified way, which had been a long-standing theoretical problem. 
© 2008 Elsevier B.V. All rights reserved. 
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In under-doped high-T c superconductors (HTSC's), 
nonmagnetic impurities (such as Zn) causes nontrivial 
widespread change of the electronic states. For example, 
both the local and the staggered spin susceptibilities are 
strongly enhanced around the impurity, within the ra- 
dius of the AF correlation length £af- This behavior was 
experimentally observed in the site-selective 89 Y NMR 
measurements for Zn-doped YBCO [1]. Moreover, a small 
concentration of Zn causes a huge residual resistivity, much 
beyond the s-wave unitary scattering limit [2] . These non- 
trivial impurity effects had been frequently considered as 
the evidence of the breakdown of the Fermi liquid state in 
under-doped HTSC's. 

Pure HTSC's samples also show various intrinsic non- 
Fermi liquid (NFL) behaviors free from impurity effects. 
Many of them had been explained based on the spin fluc- 
tuation theories like the SCR theory and the fluctuation- 
exchange (FLEX) approximation [3]. For example, Curie 
like behavior of the Hall coefficient i?n is naturally under- 
stood if one take the current vertex correction (CVC) due 
to AF fluctuations into account [4]. Anomalous transport 
phenomena in the pseudo-gap region, such as a drastic in- 
crement of the Nernst coefficient, are also reproduced by 
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taking account of CVC on the basis of the FLEX+T-matrix 
approximation [5]. 

In the present paper, we study a single impurity problem 
in a (N x TV) square lattice Hubbard model (with N = 64) . 
In general, it is not easy to obtain an appropriate solution 
for this model because two different kinds of strong inter- 
actions have to be taken into account on the same footing. 
To overcome this difficulty, we develop the Gl^-method, 
which is a powerful method of calculating the electronic 
states in real space. Based on the Gl^-method, we suc- 
ceeds in explaining nontrivial impurity effects in HTSC's 
in a unified way within the scheme of the spin fluctuation 
theory, without assuming any exotic mechanisms. 

In the present model, the self-energy S(r — r', e) is not 
a function of r — r' because of the absence of translational 
invariance. Therefore, the self-energy in the site represen- 
tation is expressed in a matrix form A^ 2 x A^ 2 , 11(e) . In the 
present model, the self-energy is divided into three parts; 
£ (e), impurity potential /, and <5S(e). The first term cor- 
responds to the self-energies only by the Hubbard interac- 
tion U , and the last term is the cross-term between U and 
I. 

In the Gl^-method [6] , S° is given by the FLEX approx- 
imation. SE itj (e) is given by T J2„ Gi,j (e + u)V/ tj (u) - S°, 
where G = ({G 00 }" 1 - S)" 1 and G 00 (e) is the noninteract- 
ing Green function. Moreover, V 1 — V" FLEX [G 7 ], which is 
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Fig. 1. Local spin susceptibility obtained by the GV^-method around 
the impurity site (at (0,0)). 

the effective interaction given by the FLEX-type diagrams 
composed of G 1 = ({G 00 }" 1 - I - E )" 1 . That is, G 1 = 
G|«5s=o- In the GV^-method, G is solved self-consistently 
whereas V 1 is a partially self-consistent function. As proved 
in ref. [6], GV^-method is much inferior to a fully self- 
consistent GU-method. 

Hereafter, we take the unitary limit / = oo. Numerical 
results of the spin susceptibility given by the GF / -method 
in site representation, x /s ( r , r )j are explained in Ref. [6]. 
Here, we present the local susceptibility x /s (r, r) in Fig. 1, 
around the impurity site (at (0,0)) along the x-direction. 
Here, we put (t,t',t") = (-1, 1/6, -1/5), where t' , t', and 
t" are the nearest, the next nearest, and the third nearest 
neighbor hopping integrals, respectively. We assume U = 8 
for hole-doped system (YBCO) and U — 5.5 for electron- 
doped system (NCCO). T = 0.02 corresponds to 80K in 
real systems. At r = (7,0), x /s ( r J r ) takes almost an bulk 
value without impurity. However, it increases gradually as 
one approach the impurity. The temperature dependence 
of x /s (r, r) away from the impurity is moderate. On the 
other hand, it increases drastically around the impurity as 
temperature decreases. At the same time, the radius of the 
enlargement of x /s ( r , r ) increases. The obtained result is 
highly consistent with NMR measurements [1] . 

Here, we discuss the physical reason for the obtained im- 
purity effect. In the FLEX approximation, the AF-order 
(in the RPA) is suppressed by thermal and quantum fluc- 
tuations. In the Gy z -method, the reduction of fluctuations 
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Fig. 2. Temperature dependence of p for NCCO with dilute impurities 
given by the GV 7 -method. 

due to an impurity gives rise to the enhancement of sus- 
ceptibility. However, this mechanism is absent in the RPA. 
Therefore, the enhancement of susceptibility is tiny within 
the RPA. 

Next, we discuss the transport phenomena in the pres- 
ence of dilute impurities on the basis of the Gl^-method 
[6] . Figure 2 shows the resistivity p for NCCO with 7ii mp = 
0, 0.01 and 0.02. The obtained result shows a huge parallel 
shift of resistivity at finite temperatures (Ap) due to im- 
purities, far beyond the s-wave unitary scattering limit. As 
T decreases, nonmagnetic impurities cause a "Kondo-like 
upturn" of p below T x , reflecting an extremely short quasi- 
particle lifetime around the impurities. We see that T x de- 
creases as riimp does. The obtained Kondo-like upturn of p 
strongly suggests that the insulating behavior of p observed 
in under-doped HTSC's is caused by residual disorder in 
the Cu02-plane, or residual apical oxygen in NCCO. 

The present study reveals that a single impurity strongly 
affects the electronic states in a wide area around the im- 
purity in the vicinity of the AF-QCP. Here, we developed 
the GT^-FLEX method, which is a powerful method to 
study the impurity effect in strongly correlated systems. 
Using the GV 1 method, characteristic impurity effects 
in under-doped HTSC's are well explained in a unified 
way, without introducing any exotic mechanisms which 
assume the breakdown of the Fermi liquid state. Qualita- 
tively, these obtained numerical results are very similar 
for YBCO, LSCO and NCCO. We succeeds in explaining 
nontrivial impurity effects in HTSC's in a unified way in 
terms of a spin fluctuation theory, which strongly suggests 
that the ground state of HTSC is a Fermi liquid. We expect 
that novel impurity effects in other metals near AF-QCP, 
such as heavy fermion systems and organic metals, could 
be explained by the GU 7 -method. 
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